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Abstract 

Dynamics of entanglement due to intensity-dependent interaction 
between a two-level atom and a single-mode electromagnetic field in a 
Kerr medium is studied. The form of the interaction is such that the 
Hamiltonian evolution is exactly solvable. The Hamiltonian is shown 
to be a deformed Jaynes-Cummings model admitting a closed, sym- 
metric algebra. Dynamics of population inversion and atom-field en- 
tanglement are studied taking the initial state of the field to be either a 
coherent state or a squeezed vacuum. Analysis is extended to the case 
of a two-mode cavity field interacting with a two-level atom. For the 
two-mode case, the initial field is a pair coherent state or a two-mode 
squeezed vacuum. Effects due to nonlinearity, intensity-dependent in- 
teraction and detuning on the dynamics are discussed and compared 
with those of the single-mode case. 

PACS: 42.50.Ct, 42.50.Dv, 03.67.Mn 
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I. Introduction 



Quantum theory allows for correlations that are not possible in any clas- 
sical stochastic process. The source of such correlations is entanglement 
PQ. Any system with two or more degrees of freedom has the possibility of 
being in an entangled state. Entangled states of atom and electromagnetic 
field in a cavity can be generated in the interaction of a two-level system 
with an external field. A well studied model of atom-field interaction is the 
Jaynes-Cummings (JC) model [21 [3j HI [5] , wherein the atom is treated as 
an electric-dipole interacting with a quantized electromagnetic field. The 
Hamiltonian is 

Hjc = ojo) a + ^i/(T^ + \{a) <7_ + aa + ), (1) 

where A is the coupling constant, uj is the field frequency and v is the 
atomic transition frequency. The operators a and a) are the annihilation 
and creation operators of the field; a z is the Pauli matrix, <t + and cr_ are 
the raising and lowering operators of the atomic states. The ground and the 
excited states of the atom are represented by \g) and |e) respectively. The 
action of the relevant operators on these states are 

a z \e) = |e), a z \g) = -\g), 
o+b) = |e), <7+|e) = 0, 
o"-|e) = \g), a--\g) = 0. 

The first two terms in the Hamiltonian Hjc correspond to the energies of 
the field and the atom respectively. The last term accounts for the atom-field 
interaction. This model has been used extensively in studying various quan- 
tum optical systems. A generalization of this model is useful in analysing 
position-dependent interaction strength [61 [8]. To account for such a 
dependence, the interaction term is modified to 

Hm = \(f(a Jf a)aa + + adjoint). (2) 

The interaction term in the JC model is obtained when f(a^a) = I, the 
identity operator. With any other form for f(a^a), the model is referred 
as Nonlinear Jaynes-Cummings(NJC) model. The operator-valued function 
/ carries information about the position-dependent interaction strength. If 
the field has nonlinear dependence on the field amplitude, as in the case of 
Kerr interaction [9l \TU[ [TT] , the Hamiltonian is generalized further to 

H = u> a) a + + X^ 2 o 2 + X(f(a^a)aaj r + adjoint). (3) 

This Hamiltonian, with f(a'a) = I, models the interaction between a field 
in a Kerr medium and a two-level atom. The quadratic term accounts 



2 



for the Kerr effect. We consider a special form of f(a f a), so chosen that 
the model Hamiltonian in Eq. (3) is exactly solvable |12[ 113] , Choosing 
f{a)a) = y/l + kd^a, where < k < 1, the Hamiltonian H is written as 

H s = ujk ] K + ^va z + \(K j fr- + adjoint), (4) 

where x = kio. The suffix "S" indicates that the Hamiltonian describes the 
interaction of a single-mode of a cavity field with a two-level atom. The 

operators k = \1 + ka^aa and = a) vT-j- ka^a are deformed annihi- 
lation and creation operators respectively. In this work, the parameter k is 

taken to be nonnegative. Nevertheless, when k = —1, the operators k and 

are the well-known Holstein-Primakoff realization of spin operators. It is 
interesting to note that when k = 0, the Hamiltonian is the JC Hamiltonian. 
Nonzero values of k amount to including Kerr effect and intensity-dependent 
interactions. In essence, nonlinearity and intensity-dependent interactions 

are included via deformation. The operators k, k' and ko = ka'a + \ 
satisfy the commutation relations, 

[k,k^] = 2k , [k ,k^} = kk\ [k ,k] = -kk, (5) 

with [ , ] denoting the commutator of the enclosed operators. Hence, the 

symmetric set {k, k\ k$} forms a closed algebra under commutation. Fur- 
ther, when k = the algebra is the Heisenberg- Weyl algebra generated by 

{d,a^,l}. Another interesting limit corresponds to k = 1. In this case, the 
algebra defined in Eq. (5) is SU(1,1) algebra. Thus, two important algebras, 
namely, the Heisenberg- Weyl and the SU(1,1), are the special cases of the 
generalized algebra defined in Eq. (5). Mathematically, choosing nonzero k 
values amounts to deforming the Heisenberg- Weyl algebra. Thus, the Hamil- 
tonian Hs can be viewed as a deformed JC Hamiltonian. Though there is 
a formal resemblance between Hjc and H$, the deformation (k ^ 0) allows 
to include nonlinearity and intensity-dependent coupling. Under suitable 
limits, the Hamiltonian Hs approximates many well-studied Hamiltonians. 
For instance, when k = 0, Hs becomes the usual JC model; when the mean 
photon number of the field satisfies k{cL i a) >> 1, then Vl + k'a^a ~ Vka^a 
and Hs approximates the Buck-Sukumar Hamiltonian [6]. 

In the single-mode field-atom bipartite system, the field is a continu- 
ous variable (infinite-dimensional Hilbert space) system while the atom is 
described in a two-dimensional Hilbert space. Recently, much work has 
been done on systems with two continuous variable subsystems; the entan- 
gled coherent states [HI [15] being the most widely studied. Such entangled 
continuous systems have desirable features which can be exploited in tele- 
portation, quantum communication, etc [TBI [M [T71 [TBI fT9l 1201 [T5l I2T1 [22] . 
Two-mode electromagnetic fields can be entangled by allowing them to in- 
teract with a two-level atom in a cavity. A basic scheme to perform this was 



3 



explored in [23] by employing a two-mode JC model. Various versions of 
the two-mode JC model [23] have been used in important contexts like the 
phase damping [25] . entanglement distribution [26], entangling field modes 
with a moving atom [27], entropy squeezing |28j etc. In this work, a natural 
extension of Hs is studied. One possibility is to consider 

H T = uxklkx + uj 2 k\k 2 + ~a z + \{k\kl&- + K 1 K 2 a + ), (6) 



where K\ = \/l + kidJa a and K 2 = y 1 + k 2 tfbb. The creation operators 

of the two modes are a and b respectively; the corresponding annihilation 
operators are and a' . The suffix "T" indicates that the composite system 
involves a two-mode field. The operators k\ and k 2 correspond to the two 
field modes of the cavity. The deformation parameters k\ and k 2 satisfy 

< ki, k 2 < 1. The operators in each of the sets {K\ , k\ , K\q = k\a)(x + ^} 
and {K 2 , k\, K 2 q = k 2 b^b+^}, satisfy the closed algebra; defined in Eq. (5). 
The dynamics dictated by the two-mode Hamiltonian Ht is exactly solv- 
able. This Hamiltonian, like its one- mode counterpart Hs, approximates, 
many of the known two-mode JC models. For instance, if k\ = k 2 = 

then Ht corr esponds to the Hamiltonian studied in the context of phase 
damping [25J. When the mean photon number is large enough so that 

1 + k{a)a) ~ k{a)a) and 1 + k(b'b) ~ k(Wb), then the interaction term 

is approximated by Vb^b V a^a aba + + adjoint and the resulting Hamilto- 
nian is the two-mode Buck-Sukumar Hamiltonian. 

In the present work, dynamics dictated by the Hamiltonian H$ and Ht 
are studied; in particular, the dynamics of population inversion and of entan- 
glement are analyzed. This paper is organized in the following way. Results 
on the single-mode system are presented first and those on the two-mode sys- 
tem are presented subsequently. In section II, dynamics of entanglement in 
the single mode case is discussed. The effects of nonlinearity and intensity- 
dependent coupling on the dynamics are studied. In section III, analyses on 
two-mode case are presented. Dynamics of various tangles in the system are 
studied, and wherever possible, the results are compared with those of the 
single mode case. 

II. Single-mode NJC model 

In the single-mode NJC model, there are two degrees of freedom in the 
system, corresponding to the field and the two-level atom respectively. The 

Hamiltonian of the system is taken to be Hs- The initial state of the atom- 
field system is taken to be a product state. Entanglement between the atom 
and the field is generated by the interaction. The time evolution of the 
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atom-field state is obtained by solving 

i^lV(*)>=^(t)>, (7) 

where \ifi(t)) represents the state of the atom-field system at time 't'. The 
explicit form of \ip(t)) is J2^=o [Ce,n(*)l e > n ) + C giTl (t)\g, n}]; the state |e,ra) 
[respectively, \g,n) ] represents the field in the number state |n) and the 
atom in the state |e) [respectively, \g)]. The coefficients C ejn (t) and C g ^ n {t) 
are [13] 



e~^C e , n {t) = [cos(^)-^sin(^)]C e ,„(0)-^ sin(^) C„ n+1 (0), 

(8) 

e iA ^C g , n+1 (t) = [cos(%^)+^sin(%^)]C s , n+1 (0)-^ sin(%^) C e>n (0), 

(9) 

where A = v— u>, A n = A—2kun, ft n = + 4A 2 ??2 and ry n = ^/(l + n)(l + kn). 
These analytical results have been used to arrive at the results presented sub- 
sequently. 

The initial state of the field is taken to be either a coherent state, 

\ a ) = e 2 E ^=l n )> a£ C, (10) 



n=0 Vn! 



or a squeezed vacuum 



\z) = -j== E(" 1 )'^r(^ tanhr)'|2/>, z = re i0 , r = \z\. (11) 



The states considered in this work correspond to a = y30 and z = 2.402. 
The values are chosen so that the mean number of photons in the CS and 
that in the SV are equal to 30. The initial state of the system is (J2 n fn\ e : n ))i 
where the coefficients f n of the field state are chosen to be those of \a) or 
\z). Thus, the initial conditions, C e>n (0) and C 9in (0), are 

C e>n (0) = f n ; C 9jn (0) = 0. (12) 

These are the initial conditions when the atom evolves from its excited state. 
The density operator for the system is the projection operator \tp(t))(ip(t)\, 
expressed in the {|e, n), \g,n)} basis as 

oo 

p = \4,(t))(m\= E E C a)n (t)Cl m (t)\a,n)(b,m\. (13) 

a,b=e,g n,m=0 
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Population inversion Ws(t) is denned as the difference in the probabili- 
ties of the atom to be in the states |e) and \g). Hence, 



oc 



W s (t) = Tr[pa z ] = ]T[|C e ,„(i)| 2 - \C g , n (t)\ 2 ] 



(14) 



n=0 



where Tr stands for tracing over both the atomic and field states. Evolution 
of Ws(t) is shown in Fig. [TJ when the field is initially the coherent state \a). 
The figures stacked in the first column correspond to k = 0; the second and 
third columns correspond to k = 10 and k = 10~ 3 respectively. Different 
rows of figures give information on the effect of detuning. The figures in the 

first row correspond to A = 0. A special feature of the Hamiltonian Hs is 
the existence of a minimum for the Rabi frequency Q n , for some n = n. The 
numerical value of n is related to the other parameters through [13] 



If the mean value of (a) a) is chosen to be n, the corresponding A satis- 
fying the Eq. (15) is the critical detuning A c . In the resonant JC model 
(k = 0,A = 0), the perfect collapses of the population inversion are Gaus- 
sian modulated cosines [29|. [30] . The dynamics of Ws (t) is known to have 
distinct features, such as 'super structures' in the "collapse - revival" pat- 
tern, when A ~ A c , i.e., the envelope of collapse or revival is not a Gaussian 
modulated cosine [9] . This is shown in the figure located at the third row of 
the second column in Fig. [TJ It is to be noted that a minimum in the Rabi 
frequency is possible only when k ^ 0. 

The first column of figures in Fig. [T] shows the effect of higher detuning 
on the dynamics of Ws (t) ■ Larger detuning inhibits the interaction between 
the atom and the field. This is expected based on the Fermi's Golden rule 
which states that interaction probability becomes smaller as the field fre- 
quency deviates from the resonant condition. Hence, with large detuning 
the atom evolves such that it has substantial overlap with its initial state. 
The last figure corresponds to A = 0.016 and the Ws(t) profile indicates 
that the atom practically remains in the excited state. 

Dynamics dictated by Hs differs from that of Hjc due to the Kerr 

term and the intensity-dependent interaction in Hs- The strength of the 

nonlinearity(oc ka) 2 a 2 ) and that of the coupling(oc VT+ ka^a) increase with 
k. It is found that if k = 0.01, even under resonant condition, the state of 
the atom has substantial overlap with the excited state |e). This implies 
that the presence of cavity field has not resulted in any significant change of 
atomic state. As noted earlier, the atom-field interaction becomes insignifi- 
cant if the interaction is non-resonant due to large detuning in the usual JC 



n = 



kA- \ 2 m 2 (l + k) 
2kw(k + A 2 ) 



(15) 
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model, which is the undeformed H$- Thus, the effect of large "fc" mimics the 
effect of large detuning. The origin of this "inhibition of emission" from the 
excited state is mostly due to the the nonlinear term ka^ 2 a 2 . In the case of 
coherent states, ka^ 2 a 2 term contributes |a| 4 + \a\ 2 to the total energy while 
the contribution from the interaction term \/T + ka^aa is ~ \/l + k\a\ 2 a. 

For small value of mean photon number ((a^a) ~ 1), the contributions of 
the nonlinear term and the interaction term to the energy are nearly equal. 
However, for large amplitude coherent state, the nonlinear term dominates 
over the interaction. To bring out these features the variation of Ws{t) is 
plotted in Fig. [3] for \a\ 2 = 1 and \a\ 2 = 30. The interaction is assumed 
to be resonant. From the figures it is seen that as k increases, the quantity 
Ws(t) stays close to unity if \a\ 2 = 30 and oscillates if \a\ 2 = 1. In the 
large amplitude case, the population inversion remains close to unity while 
it swings between 1 and 0.2 when \a\ 2 ~ 1. The former corresponds to dom- 
inant nonlinearity while the later corresponds to nearly equal magnitudes 
of the nonlinearity and the interaction. Thus, the "inhibition of emission" , 
despite the presence of a resonant field, is due to the nonlinear term. 

Even at smaller values of k the effect of nonlinearity can be dominant. 
The last column of Fig. [1] corresponds to k = 10~ 3 , an order of magnitude 
smaller than 10 -2 . The initially excited atom is left almost undisturbed in 
the resonant evident from the nearly time-independent Ws(t). As 

the detuning increases, the atom begins to interact with the field, and, con- 
sequently, the oscillations in Ws(t) are larger. This is inferred by comparing 
the figures in the first row (zero detuning) of Fig. [1] with those in the suc- 
cessive rows (higher detunings). These counter-intuitive features, that is, 
noninteracting, resonant field and interacting, non-resonant field, are due to 
the presence of nonlinearity. 

In Fig. [21 the dynamics of Ws(t) is shown when the field is initially 
a squeezed vacuum. Some of the features present in the interaction of the 
atom with coherent state are absent. Even though, the mean number of 
photons in the squeezed vacuum is chosen to satisfy Eq. (15), there are no 
"super" structures in the population inversion. However, as in the coherent 
state case, the field-atom interaction is strengthened at nonzero detuning 
when the value of the deforming parameter k is nonzero but small ( 10~ 3 ). 
The mean number of photons in the CS and the SV are the same. Hence, 
the differences in the dynamics are due to the different photon statistics of 
the two states. 

Interaction between the atom and the field entangles them. The initial 
state X)/«l e ) n ) is a product state and hence the entanglement is zero. To 
study the effectiveness of the interaction in entangling the atom and field 
system, linear entropy is used as the measure of entanglement. If p is the 
density operator of the composite system, its linear entropy (L) is defined 
as 

L = 2(1 — Tr[p 2 ]) = 2(1 — Tr[p 2 ]). (16) 
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Here pf [respectively, p a ] is the reduced density operator of the field [re- 
spectively, atom] obtained by partially tracing over the atomic [respectively, 
field] states. For instance, the reduced density operator p a for the atom is 

J2n Ce,n{t)C*^ n (t) J2n\Cg,n(t)\ J 

Even though other measures such as von Neumann entropy exist, in the case 
of pure states all entanglement measures are equivalent. In this work, linear 
entropy is used as the entanglement measure as it is easier to compute. The 
two quantities, namely, Ws and L, are related by 

L=l-Wl-4\J2 C *e,nC 9 ,n\ 2 . (18) 

The quantity | J2 Cl n C g ^ n \ is the atomic coherence since this particular term 
appears as off-diagonal element in the reduced density matrix of the atom. 
The time-evolution of population inversion shows a striking resemblance to 
the entanglement if the initial field state is a squeezed vacuum, see Figs. [2] 
andO This feature is not to be seen in the case of coherent states, see Figs. [T] 
and[5J This observation is explained by the relationship, Eqn. (18), between 
the linear entropy and the population inversion. It turns out, as indicated 
by numerical calculations, that the atomic coherence vanishes during evo- 
lution if the initial state is a squeezed vacuum. Consequently, L = 1 — Wj 
and there is a good resemblance of the L(t) and Ws(t) profiles. In the case 
of CS, the atomic coherence is not small. As a result, the L(t) and Ws(t) 
profiles do not match. Interestingly, in this case the atomic coherence spikes 
to larger values during the course of evolution. Whenever, the coherence is 
large, the reduced density matrix is less mixed. This, in turn, means that the 
entanglement of the composite system is less. This is the origin of low en- 
tanglement spikes in the case of atom-coherent state interaction, refer Fig. [H 

The time-evolution of linear entropy (L) is shown in Fig. HI when the 
cavity field is initially a coherent state. In Fig. \5\ results on the dynamics 
of L are given for the squeezed vacuum case. The first row of the figures 
correspond to A = (resonant case). The figures in the first column corre- 
spond to k = while those in the second column correspond to k = 10 -3 . 
In the resonant case, the amplitude of variation of L smaller if the value of 
k is nonzero. The effect of deformation is to minimize the effectiveness of 
the interaction in entangling the atom and the resonant field. This readily 
corroborates with the fact that the population inversion is nearly unity for 
this case, refer Fig. [T] and Higher values of population inversion (~ 1) 
imply that the atomic state is ~ |e); consequently, the composite system is 
a product state and the linear entropy (equivalently, the entanglement) is 
nearly zero. Hence, with the increase of k, the entanglement in the atom- 
field system becomes zero, even in the resonant case. In the non-resonant 
NJC model, higher values of k make the time-evolved atomic state to have 
substantial overlap with the ground state. This is in contrast with the non- 
resonant JC model in which the interaction becomes less effective in causing 
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Table 1: Mean entanglement in single- mode NJC model 
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k 


Mean entanglement 


Coherent state 


Squeezed vacuum 












0.84 


0.97 





1 


x 10" 


-4 


0.89 


0.94 





1 


x 10" 


-3 


0.04 


0.34 


0.01 









0.56 


0.64 


0.01 


1 


x 10" 


-4 


0.91 


0.77 


0.01 


1 


x 10" 


-3 


0.75 


0.26 


0.0161 









0.3 


0.44 


0.0161 


1 


x 10" 


-4 


0.74 


0.62 


0.0161 


1 


x 10" 


-3 


0.81 


0.24 



transitions as detuning increases. Thus, in the non-resonant case, the NJC 
interaction is more effective in entangling the atom and the field than the 
JC interaction. 

A simpler way of assessing the effectiveness of the interaction in entan- 
gling the field and the atom is to study the mean entanglement L, defined 
as [31] 

L = l J\{t')dt\ (19) 

is used. The value of T is chosen so that gT is 100. The values of the mean 
entanglement is shown in Table 1 for the single mode case. Small nonzero 
values of k, of the order of 10 -4 leads to better entanglement than zero or 
large values for k. In the resonant case L approaches zero as k increases 
from 10 -4 to 10 -3 . In the non-resonant case the entanglement does not 
become small even if k increases to 10 -3 . It is interesting to note that in the 
single-mode case, entanglement is more effective in the non-resonant case if 
the initial field is a CS in comparison to the case of initial field being in the 
SV. 

III. Two mode NJC model 

In this section, the interaction between a two-level atom and a two-mode 
cavity field is studied. Two types of fields are considered: pair coherent (PC) 
state |C) [32] and two- mode squeezed vacuum (TSV) [33]. The respective 
Fock state representations are: 

io = ( 2 °) 



where Nq = W'^o(2|CI); ^o(2|CI) * s the modified Bessel function of order 
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zero and ( £ C is the amplitude of the state; 



M = ^i-M 2 K n >> ( 21 ) 

n 

where /u G C is the amplitude of state and < < 1. Both these states 
involve superposition of paired states \n, n), which have equal number of 
quanta in both modes. The unpaired states \n,m)(n ^ m) are not present. 
In the quantum entanglement context, the TSV is an important example of 
entangled Gaussian state and the PC state is entangled but non-Gaussian. 
Interestingly, TSV ari ses also in the context of teleportation of uniformly 
accelerated objects [34"]. 

The Hamiltonian of the system is taken to be Ht, the two-mode exten- 
sion of Hs- The state of the system at time "t" is \^(t)) = J2^=o (Ce,n,n(i)|e, n, 
C g n ,n(t)\g,n,n)). Unpaired states are absent in the evolved state (^(i)) as 
well. This feature is due to the type of interaction being considered and 
the absence of unpaired states in the initial state. In what follows, we set 
k\ = &2 = k. With this choice for the two deformation parameters, the 

interaction term Hj = X^J (1 + A;dta)(l + ktfb)aba + + adjoint gives, 

Hj\e,n,n) = A(l + kn + k)\g, n + 1, n + 1), (22) 

Hi\g,n,n) = A(l + kn — k)\e, n — 1, n — 1). (23) 

That is, a paired state is transformed to another paired state. More precisely, 
under the action of Ht, the paired state |e,n, n) becomes the paired state 
\g, n + 1, n + 1) and vice- versa. Hence, the evolved state is a superposition of 
paired states if the initial state is a superposition of paired states. Another 
consequence is that there are invariant subspaces in the Hilbert space of the 
tripartite system. It is easily recognized from the transformations indicated 
in Eqs. (21) and (22), that the span of {|e, n, n), \g, n + 1, n + 1)}, for a fixed 

n, is invariant under the action of Hj. Moreover, the invariant subspaces are 
disjoint. These features imply that the dynamics of the coefficients C e)n ^ n (t) 
and Cg tn+ i in+ i(t) are coupled to each other and do not depend on the dy- 
namics of other coefficients. The Schrodinger equation for the dynamics of 
the coefficients C ejnin (t) and C Sin+ i jn+ i(i) is the matrix equation, 

■1 ( C e ^ n {t) \ = 

dt \ Cfl,n+l,n+l(*) J 

Here h\\ = lo\ (n + | + kn 2 — kn) + uj2{n + | + kn 2 — kn) + 4, 

/122 = u>i(n + \ + kn 2 + kn) + uJ2(n + \ + kn 2 + kn) — y, 

and h\2 = /121 = A(l + n)(l + kn) . Here, A = v — (uj\ +0J2) is the detuning. 

This matrix equation is solved to give 

= [cos( ^ ) _^ sin( ^ )]Cerin(0) 



/in 


^12 


f C e ,n,n{t) \ 


^21 


^22 





(24) 
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: — cinf _ 



sm(^f) C g , n+1>n+1 (0), (25) 



<"n,n 



! 7 Cp )n+ i in+ i(t) = [cos(— — ) + — sm(— — )]C 9>n+ i in+ i(0) 

2iXrj nn . Q. nn t fna\ 
75 — - sm(— 5— ) C e> „,„(0), (26) 



where A n , n = A - 2fcwin - 2kuj 2 n, tt n>n = „ + 4A 2 r/2 n and rj niTl = 

(1 + n)(l + fcn). These analytical solutions are used in calculating popula- 
tion inversion and entanglement measures in the tripartite system. 

The reduced density operator pa for the atomic system is 

p A = ( |Ce,n,n(i)| 2 ^n^e,n,n{^)^g,n,n{'t) j 

\ En^e|ji,nW^ lt i,iiW J2n\^g,n,n(t)\ J 

For the two-mode case, the population inversion, denoted by Wxif), is given 
by 

W T {t) = ^[|C e , n ,„(t)| 2 - \C g , n , n (t)\ 2 }. (28) 

n 

The evolution of population inversion is given in Fig. [6] for the PC state 
and in Fig. [7J for the TSV. In both the figures, the successive columns 
correspond to k = and k = 2 x 1CP 3 respectively. And the successive 
rows correspond to different detuning parameter values (A) starting from 
zero detuning. In the resonant two-mode JC model, i.e., A = and k = 
0, the Rabi frequencies satisfy 0, njTl t = gt(l + n); these frequencies are 
commensurate and so the dynamics is periodic. As a consequence, Writ) 
shows periodic features in the resonant two-photon JC model. At the time 
instances when gt = In (I, a positive integer) the atom returns to the initial 
state, unlike in the single mode case. This periodicity is destroyed in the 
non-resonant interaction as the Rabi frequencies are incommensurate. The 
atomic evolution from the excited state is increasingly inhibited with larger 
detuning which is evident from the column 1 of Figs. [6] and [Jj It is to 
be noted that the Rabi frequency O n)n depends only on a single quantum 
number "n" , a consequence of the initial field being a superposition of paired 
state. As a function of n, the two-mode Rabi frequency exhibits a minimum 

at n = n, obtained by solving g t ' n \ n =n = 0, provided k ^ 0. The detuning 
corresponding to n = h is called critical detuning, which is given as 



a 2 3 - 1 - 

A* fi = A^ = kcoN + + k)(2 + N + kN + -N 2 k) + -N 3 k% (29) 

kuj II 

where A = 2n is the total mean photon number of the two-mode radiation 
field. At critical detuning, the evolution of Wr{t) has noticeable features. 
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The envelope of Wx(t) appears periodic and this is lost when detuning 
shifts away from A' c . This is seen by comparing the figures in the third 
row(A = A^,) with the those in the second row(A < A' c ) of Fig. [UJ Further, 
nonzero k value (refer the second column of Fig. 6) significantly influences 
the envelope pattern in the population inversion. As in the case of single 
mode NJC, the critical detuning alters the population inversion dynam- 
ics qualitatively. A difference in the evolution at critical detuning, in the 
single-mode and two- mode cases is to be pointed out. In the single-mode 
case, super structures are seen if the initial field state is a CS while such 
structures are not prominent if the initial state is a SV. But in the two-mode 
case, similar structures in the envelope of Wt{£) are observed whether the 
initial field state is a PC state or a TSV. 

The system described by Hp is a tripartite system: a two-level atom 
and two field modes. The Hilbert spaces associated with the field modes 
are of infinite dimension. There is no known measure of entanglement if 
any subsystem of a tripartite system is of infinite dimension. An assessment 
of entanglement in the system is obtained by studying bipartite entangle- 
ment in different partitions of the system. In the present case, the relevant 
partitions are: (i) two-level atom as a subsystem and the two field modes 
together form the other subsystem; (ii) two-level atom along with one of the 
field modes as one subsystem and the other field mode as the other sub- 
system; and (iii) the third possibility is to exchange the roles of the fields 
in the previous case. The entanglement measures of the different bipartite 
partitions are called tangles (T). The states considered in this work are 
pure states and hence the various tangles are chosen to be the respective 
von Neumann or linear entropies. Thus, 

T AFl ® F2 = 2[1 - Tr(p 2 A )} = 2[1 - Tr(p 2 Fl ® F2 )], (30) 

Ta®f u f 2 = -Tr[p A ® Fl log 2 (pA®Fi)] = -Tr[p F2 log 2 QS F2 )], (31) 

Ta®f 2 ^ = -Tr[p A ®p 2 log 2 (pA®F 2 )] = ~Tr[p Fl log 2 Q5 Fl )]. (32) 

The suffix (A, F± £5> F2) indicates that the atom is treated as a subsystem 
and the two field modes (F\ and F2) are treated as a single entity is another 
subsystem of the relevant partition. Similar convention is adapted for the 
other suffixes. The initial field state, a superposition of paired states is 
symmetric under the exchange of photon numbers of the two modes and 
this symmetry is preserved during evolution. Consequently, Ta® Fi>F2 and 
Ta®f 2 ,f 1 are identical. In Figs. [5] and EE only Ta®Fi,f 2 is shown. The 
calculation of Ta®f-i f 2 requires the reduced density matrix of the atom and 
the first mode, which is 

PA®F 1 = \Ce,n,n\ 2 \en)(en\ + C e ^ n C* g ^ n \en){gn\ 

n 

+c g \gn)(en\ + \Cg )n , n \ 2 \gn){gn\. (33) 

Another quantity that is of interest is the entanglement between the two 
field modes. Quantum relative entropy E Fl) p 2 is a measure of entanglement 
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in a bipartite field state involving only the paired states [35]. The field 
density matrix is 

PF U F 2 = E (rMt))m)\r) (34) 

r=e,g 

= ^2an,m\n,n)(m,m\ , (35) 

n,m 

obtained by tracing over the atomic states in the complete density matrix 
of the tripartite system. The density matrix element a nm , in terms of 
the coefficients C e ^ n ^ n and 9jrt)ri , is C e ^ n ^ n C^, m m + C g ,n,nC^ m m - The von 
Neumann entropy is 

S(pf u f 2 ) = -Tr[PFi,F2 log 2 (PFi,F 2 )], (36) 
and the quantum relative entropy is 

E F lt F 2 = ~ E l °S2( a n,n) ~ S(PFi,F 2 )- (37) 
n 



The tangles Ta,FtSiF 2 , Ta®f 1: F2 and the relative entropy E(pp lt p 2 ) are 
shown in Figs. [8]and[9]fbr the system when the initial states are PC and TSV 
respectively. As in the case of Figs. [6] and the two columns correspond 
to k = and 2 x 10 -3 respectively. The three rows of figures correspond 
to A = 0, A = 0.01(< A' c ), and A = 0.0161(= A' c ) respectively. The 
mean photon number of the two modes is taken to be 3. This ensures 
that the contribution from higher photon number states are negligible and 
numerical computations are reliable. Firstly, the time evolution profile of 
the tangle Tj^p u p 2 and the relative entropy Ep 1: p 2 do not resemble that of 
the population inversion Wx(t)- However, Ta 1 f 1 ®f 2 which is linear entropy, 
is related to Wx(t) through 

W T = 1 - Ta,Fi®F 2 - 4| E Ce,n,nCg,n,n| 2 - (38) 

In spite of this relation, the profiles of Wt and Ta,f 1 ®p 2 are not similar. This 
is in contrast with the single-mode case wherein the profiles of the entangle- 
ment resemble with Ws(t), at least for the smaller values of the parameter 
k. The reason for this dissimilarity is that in the two-mode case, the atomic 
coherence | J2 Cl nn Cg )n ^ n \ 2 is comparable to Ta,Fi®f 2 and does not become 
small either for the pair coherent state or the two-mode squeezed vacuum. 

The profiles of the tangles and the relative entropy are periodic, when 
k = and A = 0, corresponding to the resonant two-mode JC model. When 
the deformation parameter k is nonzero but small (= 2 x 10 -3 ), the period- 
icity in the evolution is destroyed. As noted earlier, this feature is present 
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in the evolution of Wt too. Another interesting aspect is the negative cor- 
relation between the time-evolutions of the tangle Ta,Fi®f 2 an d the relative 
entropy E, i.e., increase (decrease) of Ta,Fi®Fz is associated with decrease 
(increase) of E. The negative correlation may be interpreted as redistri- 
bution of entanglement, from the field-field entanglement to that between 
the atom and the two modes. This feature is retained even if k is nonzero 
implying that the nonlinearity and the intensity-dependent interaction do 
not hinder the redistribution of entanglement. 

In the non-resonant case, the atom does not interact with the field ef- 
fectively and hence the initial entanglement is not expected to vary much. 
It is essential to study the cases corresponding to A = and A / A' c . 
There is an increase of amplitude in the tangle oscillation as k increases 
to 2 x 1CT 3 from zero as if the atom-field interaction is effective as in the 
resonant case. Thus, small non-zero k leads to a "resonance" effect in the 
non-resonant situation. At slightly larger values of k ~ 1CP 2 , the quantity 
Ta,Fi®f 2 becomes insignificant indicating very less entanglement between 
the atom and the two modes. However, this fact is to be expected; for such 
values of k the population inversion Writ) is nearly unity implying that the 
atomic state is ~ |e) and hence the atom is not entangled with the fields. 
The need for larger k, in comparison to k = 10~ 3 in the single-mode case, 
to disentangle the atom and the two fields is because of the smaller photon 
number (« 3) in the two-mode case. 

As in the single mode case, it is useful to define the mean of the tangles 
and the relative entropy to characterize the effectiveness of the interaction 
in entangling the systems. Definitions similar to L are adopted for the en- 
tanglement measures Ta ) Fi®f 2 and Ta®Fi,f 2 an d the relative entropy Ep 1) F 2 - 
The mean is computed over a period T such that gT = 20. From the table 
it seen that there is not as much effect on the mean values with increase of 
k when compared with the single mode case. 



Table 2: Mean tangles and Mean relative entropy in two-mode NJC model 







Mean 1 


A,F 1 ®F 2 


Mean 1 




Mean Ep x% F 2 


A 


k 


PC 


TSV 


PC 


TSV 


PC 


TSV 








0.39 


0.50 


1.89 


2.41 


1.42 


1.86 





2 x 10 -3 


0.41 


0.7 


1.88 


2.4 


1.38 


1.65 


0.01 





0.58 


0.58 


2.1 


2.52 


1.45 


1.86 


0.01 


2 x 10~ 3 


0.61 


0.7 


1.97 


2.51 


1.27 


1.74 


0.0161 





0.40 


0.40 


2.08 


2.49 


1.60 


1.98 


0.0161 


2 x 10~ 3 


0.66 


0.59 


2.09 


2.51 


1.33 


1.84 



IV. Summary 

Nonlinearity and intensity-dependent interactions affect entanglement be- 
tween an atom and a cavity field. The time-evolution due to &s, wherein 
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both nonlinearity (Kerr type) and intensity-dependent interaction (oc + intensity) 
are included, is exactly solvable. The dynamical symmetry of the Hamilto- 
nian is a deformed Heisenberg-Weyl algebra. When the deformation van- 
ishes, the algebra is the Heisenberg-Weyl algebra. The deformation can be 
continuously varied so that the algebra is SU(1,1) when the deforming pa- 
rameter is unity. The deforming parameter has the physical significance that 
the nonlinearity and the interaction are proportional to it. The Hamiltoni- 
ans considered in this work approximate many other Hamiltonians that are 
useful in modeling various quantum optical systems. 

The deformed Hamiltonian Hg makes the Rabi frequency to attain a 
minimum value which allows for critical detuning. When the cavity field is 
chosen so that the mean photon number is n and there is critical detuning 
between the field and the atom, the dynamics exhibits "super structures" 
in the evolution of population inversion. In the case of single-mode inter- 
action, with a small deformation (k m 10 -4 ) the evolution of population 
inversion and entanglement are very different from those in the usual JC 
interaction (k = 0). One of the effects of nonzero k is that the collapses 
and revivals in the population inversion are less pronounced due to the Kerr 
nonlinearity. In reality, weak nonlinearities and intensity-dependent inter- 
actions are often unavoidable. This implies that in experiments it is really 
hard to see as many distinct "revivals and collapses" as indicated by the 
usual single-mode JC model. At resonance, the effect of higher k is to in- 
hibit the evolution of the atom from its excited state. These features are 
common to both the single-mode and the two-mode NJC models. Some 
aspects of population inversion dynamics are different for the single-mode 
and two-mode models. In the undeformed, resonant case the population 
inversion Ws(t) exhibits 'collapse-revival' structure when the initial state is 
a coherent state and these features are absent if the initial field is a squeezed 
vacuum. In the two-mode case, such structures are seen when the cavity is 
either a pair coherent state or a two-mode squeezed vacuum. In the nonreso- 
nant case, the deformed Hamiltonian generates more entanglement than the 
undeformed case. In the resonant case, small deformations tend to reduce 
the effectiveness of the interaction and hence the entanglement is reduced. 
In the two-mode case, the tangle Ta < f 1 ®f 2 an d the relative entropy Ef lt p 2 
of the two fields exhibit negative correlation in their dynamics. This is in- 
terpreted as the redistribution of entanglement from that between the fields 
to that between the atom and the fields and vice versa. 
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Figure 1: Variation of population inversion versus scaled time (At) when the initial 
state is |a)|e); the field state in the coherent state \a) (a = \/30) and the atom 
in the excited state |e). The coupling constant is A = 0.001. Successive columns 
correspond to k equal to 0, 10 -4 , 10~ 3 . In the first and second columns, top most 
plot corresponds to A = and successive plots correspond to A = 0.01 (< A c ), A = 
0.016061 (= A c ). In the third column, top most plot correspond to A = and 
successive plots correspond to A = 0.05(< A c ), A = 0.061061(= A c ). 
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Figure 2: Variation of population inversion with scaled time (At). The initial 
state of the field is taken to be squeezed vacuum state \z) (z = 2.402) and that 
of the atom is the excited state |e). Successive columns correspond to k equal to 
0, 10~ 4 , 10~ 3 . In the first and second columns, top most plot corresponds to A = 
and successive plots correspond to A = 0.01(< A c ), A = 0.016061(= A c ). In the 
third column, top most plot corresponds to A = andn successive plots correspond 
to A = 0.05(< A c ), A = 0.061061(= A c ). Other parameters are the same as in 
Fig. 1. 
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Figure 3: Temporal profile of population inversion with scaled time (Xt). The 
initial state of the field is coherent state |a) and mean number of photons in the 
mode is \a\ 2 = l(left) and \a\ 2 = 30(right). The initial state of the atom is taken 
to be excited state |e). The interaction is assumed to be resonant (A = 0). The 
coupling constant is A = 0.001 and k = 10~ 2 . 
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Figure 4: Entanglement versus scaled time (At). The initial state of the field is 
the coherent state \a) (a = \/30) and that of the atom is the excited state |e). 
Successive columns correspond to k equal to 0, 10~ 4 , 10~ 3 . In the first and second 
columns, top most plot corresponds to A = and successive plots correspond 
to A = 0.01(< A c ), A = 0.016061(= A c ). In the third column, top most plot 
corresponds to A = and successive plots correspond to A = 0.05(< A c ), A = 
0.061061(= A c ). Other parameters are the same as in Fig. 1. 
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Figure 5: Entanglement versus scaled time (At). The initial state of the field is the 
squeezed vacuum state z) (z = 2.402) and that of the atom is the excited state \e). 
Successive columns correspond to k equal to 0, 10~ 4 , 10~ 3 . In the first and second 
columns, top ost plot corresponds to S — and successive plots correspond to A = 
0.01(< A c ), A = 0.016061(= A c ). In the third column, top most plot corresponds 
to A = and successive plots correspond to A = 0.05(< A c ), A = 0.061061(= A c ). 
Other parameters are the same as in Fig. 1. 
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Figure 6: Variation of population inversion with scaled time (Xt) when the initial 
state is |C)l e ); the field in the pair coherent state \Q (( = 1.778) and the atom in the 
excited state |e)_. The coupling constant is A = 0.002, the mean number of photons 
in the field is N — 3 and the critical detuning is A' c = 0.016. Successive rows 
correspond to A = 0, A = 0.01(< A' c ), A = A' c . Successive columns correspond to 
k equal to 0, 2 x 10" 3 . 
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Figure 7: Variation of population inversion with scaled time (Xt) when the initial 
state is |/x)|e); the field state in the two-mode squeezed vacuum state (fi — 1.032) 
and the atom in the excited state |e). Successive rows correspond to A = 0, A = 
0.01(< A' c ), A = A' c . Successive columns correspond to k equal to 0, 2 x fCT 3 . 
Other parameters are the same as in Fig. 6. 
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Figure 8: Evolution of various biparitite entanglement measures versus scaled time 
(At). The initial state of the field is the pair coherent state \Q (£ = 1.778) and that 
of the atom is the excited state |e). The first column corresponds to k = and 
the second column corresponds to k = 2 x 10~ 3 . The top most row corresponds 
to A = and the successive rows correspond to A = 0.01 and A = 0.0161. In 
each figures, the quantities shown are Ta,f 1 ®f 2 (bottom), E(pf 1 ,f 2 ) (middle) and 
TA®F 1 ,F 2 (t°p)- Other parameters are the same as in Fig. 6. 
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Figure 9: Evolution of various biparitite entanglement measures with time (At). 
The initial state of the field is the two- mode squeezed vacuum state (/j, = 1.032) 
and that of the atom is the excited state \e). The first column corresponds to k = 
and the second column corresponds to k = 2 x 10~ 3 . The top most row corresponds 
to A = and the successive rows correspond to A = 0.01 and A = 0.0161. In 
each figures, the quantities shown are Ta,f 1 ®f 2 (bottom), E(pp 1 ,F 2 ) (middle) and 
TA®F 1 ,F 2 (t°p)- Other parameters are the same as in Fig. 6. 
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